Sum rules and short-range correlations in nuclear matter at finite temperature 
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The nucleon spectral function in nuclear matter fulfills an energy weighted sum rule. Comparing 
two different realistic potential, these sum rules are studied for Green's functions that are derived 
self-consistently within the T matrix approximation at finite temperature. 



The microscopic study of the single-particle properties in nuclear matter requires a rigorous treatment of the nucleon- 
nucleon (NN) correlations |l|,|2j. I n fact, the strong short-range and tensor components, which are needed in realistic 
NN interactions to fit the NN scattering data, lead to important modifications of the nuclear wave function. A clear 
CN| indication of the importance of correlations is provided by the observation that a simple Hartree-Fock calculation for 
nuclear matter at the empirical saturation density using such realistic NN interactions typically results in positive 
energies rather than the empirical value of — 16MeV per nucleon 0|. 
^"5 ■ Correlations do not only manifest themselves in the bulk properties but also modify the single-particle properties 
Q\ \ in a substantial way. Several recent calculations have shown without ambiguity how the NN correlations produce a 
partial occupation of the single particle states which would be fully occupied in a mean field description and a wide 
distribution in energy of the single particle strength. These two features have also been empirically founded in the 
analysis of the (e,e'p) nucleon knock-out reactions Q. The theoretical studies have been conducted both in finite 
nuclei Q| and also in nuclear matter 0, 0, 0| • 

An optimal tool to study the single particle properties is provided by the self-consistent Green's function technique 
(SCGF) || . This method gives direct access to the single particle spectral function, which should be self-consistently 
determined at the same time than the effective interactions between the nucleons in the medium. Enormous progress 
in the SCGF applications to nuclear matter have been reported in the last years, both at zero Q and finite temperature 

The efforts at T = have mainly been addressed to provide the appropriate theoretical background for the interpreta- 
tion of the (e,e'p) experiments while the investigation at finite T are mainly oriented to describe the nuclear medium 
in astrophysical environments or to the interpretation of the dynamics of heavy ion collisions. 

In any case, the key quantity is the single-particle spectral function, i.e. the distribution of strength in energy when 
one adds or removes a particle of the system with a given momentum. A possible way to analyze the single-particle 
L* ' spectral function is by means of the energy weighted sum rules. They are well established in the literature and have 
been numerically analyzed in the case of zero temperature |12j . 

The analysis of the energy weighted sum-rules can give useful insights not only on the numerical accuracy of the 
many-body approach used to calculate them but also can help to understand the properties and structure of the NN 
potential. 

This paper is devoted to study the physical implications of the fulfillment of these sum rules for single particle spectral 
functions in nuclear matter at finite T. This investigation is based on the framework of SCGF employing a fully self- 
consistent ladder approximation in which the complete spectral function has been used to describe the intermediate 
states in the Galistkii-Feynman equation. 

After a brief summary of the definitions of the single-particle spectral function, we give a simple derivation of the 
sum rules. Then we analyze the results for two types of realistic potentials, the CDBONN and the Argonne V18, 
and discuss the different behaviors based on the different strength of the short-range and tensor components of both 
potentials. 

For a given Hamiltonian H , the Green's function for a system at finite temperature can be defined in a grand-canonical 
formulation: 

ig(kt;k't') = Tr{pT[a k (t)al(t')]}. (1) 

T is the time ordering operator that acts on a product of Heisenberg field operators a,k(t) = e ltH a^e^ ltH in such a 
way that the field operator with the largest time argument t is put to the left. The trace is to be taken over all energy 
eigenstates and all particle number eigenstates of the many body system, weighted by the statistical operator, 

p= l (2 ) 
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j3 and fi denote the inverse temperature and the chemical potential, respectively. N is the operator that counts the 
total number of particles in the system, 

w = E°kW«kW (3) 

k 

N is independent of time, since it commutes with H. The normalization factor in Eq. (J2J is given by the grand 
partition function of statistical mechanics, 

Z = Tve-^ H -^ N K (4) 

For a homogeneous system, the Green's function is diagonal in momentum space and depends only on the absolute 
value of k and on the difference r = t' — t. Starting from the definition of the Green's function, we first focus on the 
case r > 0. In order to recover the expression for the ensemble average of the occupation number n(k) for r = + , the 
following definition of the correlation function, g < , includes an additional factor of — i with respect to the definition 
of the Green's function g, 

g<(k,r) = Tr{pe^ H ale-^ H a k }. (5) 
g < {k 1 r) can be expressed as a Fourier integral over all frequencies, 

/ + OC J 
— e-^A<{k,w) (6) 

if A < (k,uj) is defined by ^ 

A<(k,uj) = 2irJ2 ^ \{* n \a k \9 m )\ S(u-(E m -E n )). (7) 

nm 

This can be easily checked by inserting the eigenstates \^f n ) into the expression of the trace in Eq. JSJ. It is important 
to note that \^ n ) are simultaneous eigenstates of both the number operator and the Hamiltonian. 
A similar analysis can be conducted for r < 0, yielding a function 

A> (k, oj) = e ^-^A< (k, u). (8) 

The spectral function at finite temperature is defined as the sum of the two positive functions, A < and A > , 

A(k,u) = A<(k,uj) + A>(k,oj). (9) 

Expression (7J, for A< , can be compared to the result for the hole spectral function at zero temperature, that was 
reported in Ref [T^. 

A h (k,oj) = 2nJ2 IK^KI^)! 2 ^ - (i# - Et 1 )), (10) 

n 

where |^^) is the ground state of an A particle system and l^^ -1 ) labels the excited energy eigenstates of a system 
that contains one particle less. The physical interpretation of the hole spectral function in a system at zero temperature 
is the following: Ah(k,Lu) is the probability to remove a particle from the ground state of the ^4-body system, such, 
that the residual system is left with an excitation energy E^~ x = Eq — u). Eq is the ground state energy of the A 
particle system. It is clear that the lowest possible energy of the final state is the ground state energy of the A — 1 
particle system, so that there is an upper limit for the hole spectral function at uj = Eq — Eq~^ = /i. In a similar 
fashion, the particle spectral function A p can be defined as the probability to attach a further nucleon to the system in 
such a way, that the excitation energy of the compound system with respect to the ground state energy of the initial 
system is uj — E A+1 — Eq. In this case, one can argue that, to add a further particle, one has to pay at least the 
chemical potential, so that fi is a lower bound for uj. At zero temperature, this behavior causes a complete separation 
of the particle and the hole spectral function. 

The situation is quite different in a grand-canonical formulation at finite temperature To illustrate these changes, the 
full spectral function A, as well as A > and A K are shown in Fig. ^for three momenta around the Fermi momentum 
of a zero temperature system at the same density, p = 0.2 fm~ 3 . Numerical values for the integrated strength of A < 
are listed in Tabled Since thermally excited states \^ m ) are always included in the grand-canonical ensemble average 
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according to their weight factor e~^ Bm ~^ Nm \ one can take out a particle from a thermally excited state and end up 
in a weakly excited state close to the ground state of the residual system. This leads to a contribution to A < for an 
energy uo larger than /i. Also one has to keep in mind that we are considering a grand-canonical average. This implies 
that with the appropriate weight one also considers systems with a density larger than the mean value. For those 
systems, a removal of particles from states with single-particle energies above /i will also be possible. 
Similarly, a particle can be added to a thermally excited state, leaving the compound system in a state close to its 
ground state, so that A > (k,uo) extends to the region below fi. In any case, there is no longer a separation between 
A > and A < , and the maxima of both functions can even coincide. This is also quite obvious from the relation JHJ. 
For the T matrix approximation to the self energy reported in ^lj, one can determine the single- particle Green's 
function as the solution of Dyson's equation for any complex value of the frequency variable z, 



j(k, 



(ii) 



Using the analytical properties of the finite temperature Green's function along the imaginary time axis, an important 
relation between the spectral function and the Green's function can be derived and analytically continued to slightly 
complex values [l^: 



: +oc duj' A(k,u) 



2tt uo — to' + ir; 
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One can extract sum rules from the asymptotic behavior at large uo by expanding the real part of both expressions 
for the Green's function, Eqs. and Q12JI. in powers of — . This yields 



Re g(fc, uo) 
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Re g(k, uo) = - I [ duo' A(k, oo 1 ) + - [ 
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By comparing the first two expansion coefficients, one finds the tuq and the mi sum rules, 

'■ + 00 



L ^ (fc ' w)=1 ' 



and 



duo 
2^ 



A(k, uo)uo 



*1 
2m 



lim ReS(fc,w) 



(13) 



(14) 



(15) 



(16) 



Similar sum rules can be obtained from the higher order terms, as it was done in Ref. [T^ for m,2. Thinking of an 
arbitrary approximation scheme for to), it might be interesting to ask whether or to what extend such a scheme 
fulfills the sum rules. This is, however, not the point we want to address in this paper. In the T matrix approximation, 
the real part of the self energy can be computed from the imaginary part, using a dispersion relation, 



Re£(fc,w) = E c 



7T J-< 



dX 



ImS(fc, A + ir/) 
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(17) 



In the derivation of Eq. I|17|l. the spectral decomposition of the Green's function was already used, so it is a property of 
the T matrix approach that it automatically fulfills the sum rules. Nevertheless, besides providing a useful consistency 
check for the numerics, it is interesting to use the sum rules to compare the importance of short-range correlations 
for different realistic potentials on a quantitative level. The first term on the right hand side of Eq. I|17|) is the energy 
independent part of the self energy, 



E°°(fc) = 



d 3 fc' 
(2^rj3 



{kk'\V\W) A n(k'). 



(18) 



which can be identified with lim^^oo ReS(fc,o;), since the dispersive part decays like — for uo — > ±oo. Eq. (|18fl looks 
like a Hartree-Fock potential, however, n{k) is the momentum distribution that is determined from a non-trivial 
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spectral function A K in Eq. JBJ), assuming r = 0. In contrast, the Hartree-Fock self energy at finite temperature 
must be determined from an energy spectrum e(k) and a momentum distribution nHp{k) = /(e(fc)), where f(oj) 
is the Fermi function. Unlike n#i?(fc), the non-trivial n(k) accounts for depletion effects of the bound states due 
to short-range correlations. In this sense, S°° is a generalization of a Hartree-Fock potential. Fig. [21 illustrates the 
difference between the two pictures with the corresponding Feynman diagrams. 

All results in this paper have been obtained using the iteration procedure that was described in Ref. [H). Fully 
self-consistent spectral functions were calculated for two realistic potentials, the stiffer Argonne V18 and the softer 
CDBONN. 

The mo sum rule is fulfilled better than 0.1% in the complete momentum range. Results for the mi sum rule are given 
in Fig |3 for a temperature of T = lOMeV and a density of p = 0.2 fnU 3 . It is satisfied better than 1%. Both right 
hand side and left hand side are plotted, but the curves lie on top of each other and cannot be distinguished (solid 
lines). The lower dash-dotted line shows the mi-contribution from A K , which is always negative and goes to zero for 
high momenta, since there A < is strongly suppressed. The probability to remove a high-momentum particle from the 
system is simply very small. The upper dash-dotted line displays the contribution from A > . Due to the short-range 
correlations, there is a high-energy tail present in the spectral function, and so this contribution is already positive 
at low momenta, furthermore, it is nearly constant in this range, reflecting the fact that the high energy strength 
distribution is momentum independent. As soon as the quasiparticle peak of the spectral function is located at energies 
greater than /i, the A > contribution increases steadily, following the position of this peak. Both contributions add up 
to TOi. It is interesting to remind the fact that for free particles, the sum rules are automatically fulfilled. In this case, 
A < and A > are delta peaks that are located at the same position and their strength adds up to one. Their relative 

,2 

strength is given by the ratio of the phase space factors f(e(k)) and [1 — f(e(k))\, respectively, where e(fc) = j— in 
the free case. 

The results in Fig. [3] shows that the sum rule mi is rather sensitive to the differences in the NN potentials. The 
mi results for the CDBONN interaction is about 65MeV more attractive than the Argonne V18 result. A closer 
examination shows that this is predominantly due to the A > contribution, which is almost 50 MeV more repulsive for 
the Argonne V18. This means that the Argonne potential produces more correlations in the sense that the strength 
that effects mi is redistributed to higher energies. 

The dotted lines are the simple Hartree-Fock estimate of mi for the same temperature and density. For both potentials, 
the Hartree-Fock result makes up quite a good approximation to the sum rule. This result is interesting, since it permits 
a quantitative estimate of the amount of correlations produced by any given NN potential without a sophisticated 
many-body calculation. 

Fig.^Jreports the exhaustion of the sum rules mo (left panel) and mi (right panel) versus the upper integration limit 
u> for a momentum of k = 500 MeV. At this momentum, the quasiparticle peak is located around 100 MeV. For both 
interactions that were considered, the main contribution to mo, more than 80%, come from the quasiparticle peak 
of the spectral function. In the region far above the peak, the CDBONN saturates considerably faster. In Table. [HJ 
the upper integration limits that have to be chosen to exhaust the sum rule to a given percentage are reported for 
mo and mi and compared for both potentials. In the case of mo and the stiffer Argonne V18, one must integrate 
almost twice as far as for the softer CDBONN. The saturation of the mi sum rule is different, because in this case, 
a somewhat higher energy region of the spectral function is probed. In the right panel of Fig. one can observe 
that the quasiparticle peak contributes less than 50% to mi, and the high-energy tail becomes much more important, 
since it is weighted by a factor of w. While both potentials behave qualitatively similar up to an integration limit 
of about 700 or 800 MeV, where mi is already exhausted by about 75% for the CDBONN potential (cf.Tab. HTJ), a 
large contribution of about 40% is still above this energy in the case of the Argonne V18. One can also note, that, 
to exhaust mi completely, one has to integrate up to higher energies in the case of the CDBONN. However, these 
contributions to the spectral function above u) « 4000 MeV are weak and yield no further repulsion. 
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k [MeV] below fi [%] above fi [%} n(k) 



230 


98 


2 


0.706 


275 


77 


23 


0.481 


320 


33 


67 


0.191 


400 


71 


29 


0.025 


500 


95 


5 


0.006 



TABLE I: Strength distribution of A < . The numbers give the fraction of the integrated strength above and below the chemical 
potential /i. The last column reports the occupation number of the respective state. The parameters are the same as in Fig.Q 



% saturation m CDB m V18 mi CDB mi V18 



60 






277 


690 


75 






790 


1518 


90 


215 


311 


2250 


2860 


95 


403 


725 


3756 


3740 


99 


1388 


2277 


8545 


5720 



TABLE II: Upper integration limits of the running integrals that must be chosen to exhaust the sum rule mo and mi up to 
the fraction given in the first column. The parameters are the same as in Fig. 2] 



[6] A. Ramos, A. Polls and W. H. Dickhoff, Nucl. Phys. A 503, 1 (1989). 

[7] Y. Dewulf, W. H. Dickhoff, D. Van Neck, E. R. Stoddard and M. Waroquier, Phys. Rev. Let. 90, 152501 (2003). 

[8] W. H. Dickhoff and H. Miither, Rep. Prog. Phys. 55, 1947 (1992). 

[9] P. Bozek, Phys. Rev. C 59, 2619 (1999). 
[10] P. Bozek, Phys. Rev. C 65, 054306 (2002). 
[11] T. Frick and H. Miither, Phys. Rev. C 68, 034310 (2003). 

[12] A. Polls, A. Ramos, J. Ventura, S. Amari and W. H. Dickhoff, Phys. Rev. C 49, 3050 (1994). 

[13] W. D. Kraeft, D. Kremp, W. Ebeling and G. Ropke, Quantum Statistics of Charged Particle Systems (Akademie-Verlag, 
Berlin, 1986). 

[14] L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, New York, 1962). 



0.08 r 
0.07 
0.06 
"a 0.05 ■ 
? 0.04- 
< 0.03 
0.02 
0.01 ■ 
■ 



k = 230 MeV 




-50 V- 



k = 275 MeV 



k = 320 MeV 



-50 V- 
co [MeV] 



-50 V- 



— I 0.08 

- 0.07 

- 0.06 

- 0.05 

- 0.04 

- 0.03 

- 0.02 
0.01 

: 



FIG. 1: Spectral function for a density of p — 0.2 fm 3 and a temperature of T = 10 MeV (solid line). Various momenta are 
considered as indicated in the three panels. A < (dashed line) and A > (dash-dotted line) are also displayed. 
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FIG. 3: Illustration of the energy weighted sum rule mi (solid lines) for the CDBONN potential (left panel) and the Argonne 
V18 potential (right panel). Both right hand side and left hand side are displayed, but the sum rule is so well fulfilled that 
they are on top of each other. The contribution to mi that comes from A > and A < is indicated by the upper and the lower 
dash-dotted lines, the latter approaching zero rapidly for high momenta. The dotted line is the Hartree-Fock single particle 
spectrum. Density and temperature are the same as in Fig. 




FIG. 4: Saturation of the sum rules mo (right panel) and mi (left panel) for the CDBONN potential (solid line) and the 
Argonne V18 potential (dashed line). The momentum is k = 500 MeV. Again, temperature and density are the same as in 

Fig. [fl 



